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Using a density matrix approach to describe the process of coalescence, we calculate the coalescence 
probabilities and invariant momentum spectra for deuterons and antideuterons. We evaluate our 
expressions with a hydrodynamically motivated parametrization for the source at freeze-out which 
implements rapid collective expansion of the collision zone formed in heavy ion collisions. We find 
that the coalescence process is governed by the same lengths of homogeneity which can be extracted 
from HBT interferometry. They appear in the absolute cluster yield via an effective volume factor 
as well as in a quantum mechanical correction factor which accounts for the internal structure of 
the deuteron cluster. Our analysis provides a new interpretation for the parameters in the popular 
phenomenological coalescence model and for the effective overlap volume in Hagedorn's model for 
cluster production in pp collisions. Using source parameters extracted from a recent HBT analysis 
of two-pion correlations, we successfully describe deuteron and antideuteron production data from 
Pb-I-Pb coUisions at the CERN SPS as measured by the NA44 and NA52 collaborations. We also 
confirm the recent finding by PoUeri et al. that the different measured slopes of nucleon and deuteron 
transverse mass spectra require a transverse density profile of the source which is closer to a box 
than to a Gaussian shape. 

PACS numbers; 25.75.-q, 25.75.Ld, 25.75.Dw, 25.75.Gz 



I. INTRODUCTION 

Ultra-relativistic heavy ion collisions are used to study 
in the laboratory nuclear matter at extreme energy den- 
sity and temperature. For such conditions lattice QCD 
predicts the transition from a hadron gas to a quark- 
gluon plasma. The only observables of the hot, central 
reaction zone (fireball) created in the heavy ion collision 
are the energies and the momenta of the produced parti- 
cles. Most of these are hadrons vifhich, due to their strong 
interactions, decouple only late from the collision fireball 
and therefore carry direct information only about this 
so-called "freeze-out stage" . Direct signals from the ear- 
lier (presumably much hotter and denser) stages of the 
collision are carried by electromagnetic (photons, dilep- 
tons) and "hard" probes (jets, J/ip) which are, however, 
more difficult to study due to their much smaller cross 
sections. Developing a clear and unique picture of the 
interesting early stages of ultra-relativistic heavy ion col- 
lisions will thus require the combination of such direct 
signatures with an extrapolation backwards in time of 
the information extracted from the spectra of the bulk 
of produced hadrons. For such a back-extrapolation to 
be useful, the analysis of the measured hadron spectra 
should permit a more or less complete reconstruction of 
both the geometry and the dynamics of the source at the 
point of freeze-out. 

For this reason the last few years have seen a growing 
effort to understand the final state of heavy ion colli- 
sions quantitatively from the measured hadron spectra. 
In this project a crucial role is played by two-particle mo- 
mentum correlations. Such correlations, whether due to 
final state interactions or quantum statistics, are sensi- 



tive to the phase-space distribution of the particles in the 
source at freeze-out and thus carry information not only 
on the momentum-space structure of the source, but also 
on its size and shape in coordinate space. In particular, 
they are sensitive to correlations between the momenta 
and positions of the particles at freeze-out such as those 
which are generated by the collective expansion of the 
hot and dense collision zone. This dynamical aspect is 
very important since, once identified and quantitatively 
analyzed, it carries valuable information about the early 
equilibration processes and the resulting build-up of ther- 
modynamical (isotropic) pressure in the reaction zone, as 
well as on the time-integrated action of the equation of 
state of the quark-gluon and/or hadronic matter during 
the expansion stage. 

The production of deuterons, antideuterons and larger 
nuclear clusters via coalescence of (anti-)nucleons at 
freeze-out can be viewed as a particular type of phase- 
space correlation among the particles in the final state. 
The similarity of the physics of coalescence and of other 
types of final state momentum correlations was stressed 
before by Mrowczyhski |^ who pointed out that simi- 
lar source information can be extracted from the anal- 
ysis of both types of correlations. In the present paper 
we carry the analogy significantly farther by extending 
the investigation to systems with strong collective dyna- 
mics. We show in particular that the same "lengths of 
homogeneity in the source" Q which can be extracted 
from two-particle Bose-Einstein correlations (Hanbury 
Brown ~ Twiss (HBT) interferometry [^) determine the 
coalescence probability in the cluster formation process. 
This means that the effects of source expansion enter the 
"HBT radii" and the cluster yields in a similar way. For 
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this reason nuclear cluster spectra provide an important 
and complementary check for the source reconstruction 
from single-particle spectra and two-particle momentum 
correlations of elementary hadrons. 

As a point in case we demonstrate that we can success- 
fully reproduce the measured cluster yields and spectra 
for deuterons, antideutcrons, tritons and '^He in Pb-|-Pb 
collisions at the CERN SPS @,|-|| with a source model 
whose parameters were determined P,p^ from a com- 
bined analysis of pion spectra and two-pion correlations. 
We further show that the shape (inverse slope) of the 
measured deuteron transverse mass spectra [HJsI, when 
compared with that of the proton spectra, yields addi- 
tional constraints on the shape of the transverse source 
distribution in space which remove a specific ambiguity 
left open by the analysis of pion spectra and correlations. 

Historically, cluster formation has been characterized 
in terms of the "invariant coalescence factor" Ba, defined 
through the invariant cluster momentum spectra via the 
equation 
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It relates the cluster spectrum to the invariant momen- 
tum spectra of the coalescing nucleons at the same veloc- 
ity. The correct physical interpretation of Ba has been 
a long-standing problem ||l|,|ll|-g6| . The relative fragility 
of the nuclear clusters implies that only nucleons with 
small relative momenta contribute to cluster formation; 
this led early researchers |]lT[-p^ to an interpretation 
of Ba in terms of a momentum space coalescence vol- 
ume, parametrized in terms of a maximal relative mo- 
mentum po between the coalescing nucleons. This inter- 
pretation appeared to be confirmed by the approximate 
constancy of the Ba values observed in heavy ion colli- 
sions at the BEVALAC with beam momenta up to about 
1 GeV/nucleon, independent of the beam energy and the 
size of the colliding nuclei. In the sudden approximation 
and the thermodynamic models, on the other 
hand, Ba V^~^ is inversely related to the fireball vol- 
ume in coordinate space. In such a picture the decrease of 
the Ba values with increasing beam energy, observed in 
nuclear collisions with large ions at beam energies above 
1 GeV/nucleon (for a recent compilation see |2^), can 
be easily understood in terms of collective expansion of 
the collision zone before break-up. Good reviews of these 
early approaches can be found in |l|,|l|,|6|. 

Later work began to stress more explicitly the phase- 
space aspects of the coalescence process, starting from 
a quantum mechanical approach based on the density 
matrix of the source or the equivalent Wigner function 
formalism |l7|Jl^ , po| , p^ , ^ and its classical phase-space 
analogues [ pl| , ^4[ . In this approach the size of the cluster 
itself enters as an additional dimensionful quantity into 
the calculation of Ba, and it allows to address the prob- 



lem of energy-momentum conservation in the coalescence 
process. 

None of the existing model calculations accounts, how- 
ever, properly for dynamical expansion of the source and 
the resulting correlations between the momenta and po- 
sitions of the particles at freeze-out. These correlations 
can be included numerically by applying a density ma- 
trix or Wigner function based "coalescence afterburner" 
to the output of classical microscopic phase-space simu- 
lations like the Intranuclear Cascade [|8j, ARC or 
RQMD [ p6p9|j30[ | . While such numerical simulations may 
succeed or fail to reproduce the experimental data, in nei- 
ther case they provide a clear conceptual understanding 
of the physics entering into the calculation of Ba- In 
this paper we present an analytical approach, based on 
an explicit source parametrization, which provides such 
an understanding. Like the static model calculations be- 
fore, it is based on a quantum mechanical density matrix 
(Wigner function) approach. It thus allows for a proper 
treatment of energy-momentum conservation and to ex- 
hibit the dependence of Ba on the phase-space structure 
of the source, i.e. on its longitudinal and transverse size 
and collective flow, and on the internal cluster structure. 

The paper is organized as follows: In Sect. |l| we re- 
view shortly the classical thermal-|-flow model approach 
to cluster production in which clusters are viewed as 
pointlike elementary particles without internal structure 
which are produced by thermal emission from an expand- 
ing source. While poorly justified theoretically, this ap- 
proach has enjoyed considerable phenomenological suc- 
cess. The results of Sect, g will serve as a benchmark for 
the discussion of the quantum mechanical density matrix 
approach presented in the following sections. We will 
discover strong formal similarities which simultaneously 
provide a theoretical justification for the phenomenolog- 
ical success of the thermal approach and yield an inter- 
pretation of the thermal model parameters for the cluster 
spectra in terms of the phase-space characteristics of the 
underlying nucleon source. This is described in Sects. [II 
and [V where we discuss the quantum mechanics of the 



coalescence process and calculate a quantum mechanical 
correction factor for the classical thermal cluster spec- 
trum. We show that the latter is given by a very simple 
expression involving only the size of the cluster and the 
"homogeneity radii" of the source which can also be mea- 
sured with HBT interferometry. In Sect. ^ we extend this 
discussion from two-nucleon clusters (d a nd d) to three- 
nucleon clusters, t and ^He. In Sect. ^ we express the 
B2 value in terms of of the HBT radii and show how 
cluster yields can be used for an analysis of the chem- 
ical composition of the fireball at the point of thermal 
nucleon freeze-out. This complements the chemical anal- 
ysis of elementary hadron yields which determines the 
(usually much earlier) point of chemical freeze-out. A 
comparison of our results with heavy ion data is pre- 
sented in Sect. VII, and our conclusions are presented in 
Sect. [VH^ . 

We will use natural units h = c = = I- We de- 



2 



note by m, rrit, and pt the rest mass, transverse mass, 
and transverse momentum of a nucleon. M, Mt, and Pt 
denote the respective variables of a cluster. For a cluster 
of A nucleons we have in good approximation M = Am, 
P = Ap, Pt — Apt, and Mt — Anit. We denote by 
Y = y ~ ycra the longitudinal rapidity of a particle rela- 
tive to cm. frame of the fireball. 



II. THERMAL CLUSTER SPECTRA 

Clusters are no elementary hadrons. Their binding en- 
ergies {B = — 2.25MeV for a deuteron) are small com- 
pared to typical collision energies in the fireball created 
in a heavy ion collision, and this makes them very fragile 
objects. While the spectra of elementary hadrons de- 
coupling from a thermalized fireball along a freeze-out 
hypersurface Sf(i?) with normal 4-vector <P(j^{R) can 
be described by the Cooper- Frye formula M , 
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where fi{R, P) is a local equilibrium distribution 
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with local temperature T{R), local chemical potential 
ljLi{R), and local flow 4- velocity u^{R) (normalized to 
u ■ u = 1), a similar expression for nuclear clusters can- 
not be justified: clusters do not pre-exist as particles 
with thermalized momentum distributions in the collision 
fireball, but are only created by final state interactions 
among the nucleons during the freeze-out process. 

The absence of preformed nuclear clusters in the col- 
lision zone of a relativistic heavy ion collision at AGS 
or SPS energies is easily seen: Nuclear fragmentation of 
the target and projectile nuclei can be excluded as the 
origin of cluster production if clusters made from anti- 
nucleons (e.g. d or ■^H) are selected. For clusters with 
positive baryon number near the midrapidity region it 
was estimated in that the probability for a nuclear 
fragment to absorb the necessary longitudinal momen- 
tum transfer of several GeV to several 10 GeV per nu- 
cleon without getting destroyed can be neglected. On the 
other hand, by applying Levinson's theorem |Q to the 
two-particle terms in the virial expansion of the grand 
canonical potential, it was shown in | |3^ , |36[ | that thermal 
cluster production inside the fireball is exponentially sup- 
pressed when the binding energy of the cluster is much 
smaller than the fireball temperature. 

In a very simple-minded approach to cluster forma- 
tion at nucleon freeze-out one may postulate that a clus- 
ter oi A = Z + N nucleons with total momentum P is 
emitted whenever Z protons and N neutrons of identical 
momenta P/A happen to be at the same place R. This 
results in a generalized Cooper- Frye formula for clusters: 
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(2.3) 



where /p^n are given by Eq. ([2.21 ). This formula agrees, 
up to a quantum mechanical correction factor which will 
be discussed below, with the Hagedorn model for cluster 
production in high energy pp collisions ■ Iii spite of its 
naivity it has been quite successful [^p8[ . It will be one 
of the main points of the present paper to explain why 
this is so. For this purpose we first analyze Eq. ( |2.3|) for 
the specific source model which will later form the basis 
of our quantitative comparison with data. 



A. Analytical model for an expanding fireball 

There is mounting evidence from phenomenological 
analyses of existing data as well as from mi- 

croscopic simulations of the collision dynamics [|oj that 
during the early stages of a relativistic heavy ion colli- 
sion the reaction zone approaches a state of local ther- 
mal equilibrium. The resulting thermal pressure causes 
a collective expansion of the system via hydrodynamic 
flow. It thus makes sense to parametrize the momentum 
distr ibut ions of the hadrons just before decoupling by 
Eq. (2_^), with the local flow velocity 4-vector u{R) de- 
scribing the average particle velocity at point R and T{R) 
parametrizing the random ( "thermal" ) local momentum 
fluctuations around their average value = mu'^{R). 
The local chemical potential fii {R) and the local fugacity 
Xi{R) = exp[fj,i{R)/T{R)] parametrize the local density 
of particle species i at point R. 

For heavy particles (nucleons and nuclear clu ster s) at 
chemical potentials /i^ ^ (see below) Eq. ( |2.2| ) can 
be very well approximated by a local Boltzmann dis- 
tribution. For simplicity we will assume freeze-out at 
constant temperature, Tf{R) — T — const. (General 
arguments based on the kinetics of the freeze-out pro- 
cess show that this is generally a good approx- 
imation.) The fugacity Xi{R) is split into a constant, 
particle-specific term Ai — exp(fii/T), and a common 
density profile H{R) for all particle species. This im- 
plements the assumptions of local chemical equilibrium 
among the various particle species at freeze-out and si- 
multaneous freeze-out of all particle species which seem 
reasonable for nucleons and nuclear clusters. 

While in Pb+Pb collisions in the initial state neutrons 
outnumber protons by a factor of 1.54, the ratio n/p or 
n/p in the fireball at freeze-out is unknown. Due to parti- 
cle production an unknown fraction of the initial isospin 
asymmetry in the nucleon sector may be transferred to 
other particle species. We introduce separate chemical 
potentials fip and /in for protons and neutrons and define 
the chemical potential of a cluster by ha = Zji^ + Nfij^. 
The values of /ip and /Xn must be extracted from a fit to 
the data. 
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We consider only very central (impact parameter b « 
0) collisions. The fireball is then azimutally symmetric 
with respect to the beam axis ("longitudinal" or z axis), 
and the transverse coordinates are conveniently chosen 
as p = y/x^ + y^ and the azimuthal angle (p. Ultrarel- 
ativistic kinematics in the beam direction suggests the 
longitudinal proper time r = \Jt'^ — and the longitu- 
dinal space-time-rapidity r\ — arth(z/t) as appropriate 
longitudinal and temporal coordinates: 



= [Tchrj, p COS0, p sm(f>, rshr;) 



(2.4) 



The particle momenta are parametrized by their rapid- 
ity Y = aith {Pz/E) along the beam direction and their 
transverse mass Mt 



Pi- 



p/^ = (Mtchy, Pt cos$, Pt sin$, MtshY) 



(2.5) 



The flow 4- velocity is conveniently parametrized in terms 
of longitudinal and transverse flow rapidities rji and 774, 
respectively: 

u^{R) ~ [chrjidn-qn shr/t cos0, shr/t sini/), shyy; ch-qt) , 

(2.6) 

where thrji = Vi, i = l,t, defines the corresponding flow 
velocities. In the spirit of Bjorken we assume a scal- 
ing velocity profile vi — z/t in the beam direction while 
taking a power-law rapidity profile in the transverse di- 
rection which is independent of z and t: 



m{T,r],p) = 77, 



VtiT,ri,p) = Vf [^j 



(2.7) 



Here Ap characterizes the transverse size of the fireball 
(see below), while 77/ represents the strength of the trans- 
verse flow; the power a of the transverse flow profile is 
generally chosen as a = 1, except for some tests with 
a = 0.5 and a = 2 as noted in the text. 

As the fireball expands the scattering rate of the parti- 
cles decreases until finally the thermalization of the sys- 
tem breaks down and the particles freeze out. Consis- 
tently with the above Ansatz for the expansion flow pro- 
flle we assume that this happens at a flxed longitudinal 
proper time tq and set H{R) — H{r],p)S{T — tq). For 
the longitudinal and transverse shape of the density pro- 
flle H{R) we take Gaussians with widths A77 and Ap, 
respectively. 

W ith these ingredients the distribution functions in 
(^.3[) take the form 



/.(i?,P) = e'^'/^e-^-(«)/^i7(i?), ^ 



p,n 
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(2.8a) 
(2.8b) 



The density profile is normalized to a total covariant 
freeze-out volume Vcov- 



d^R H{R) = (2^) 5 ( Ap)2 ( A77) TO , (2.9) 



where d^i? = rdr pdpdr] d(j). For freeze-out at cons tan t 
longitudinal proper time, the integration measure in ( |2.3| ) 
over the freeze-out hypersurface is given by P - d^<j{R) — 
To Mt p dp ch (77 — Y) drj dcf). 



B. Non-zero emission duration 

In instead of a Cooper- Frye integral over a 3-dim 
freeze-out hypersurface, invariant spectra are calculated 
as a space-time integral / d'^RSi{R, P) over an emission 
function 

5,(i?, P) = Mt ch (,y - Y) 



(27r)3 



xi/(77,p) J(t). 



(2.10) 



[H{f],p) = y/(2/njH{r],p) differs from (U) only by 



the normalization. The present choice is more convenient 
for us because it absorbs some constant terms in the clus- 
ter spectra below which would otherwise scale with the 
nucleon number. However, it affects the interpretation 
of the total fireball volume {Vcov = \/ '2 / ttVcov) , and of 
the fugacity factor exp(/i/r). In the case of Ref. fj,/T 
is the fugacity averaged over the fireball; in the present 
case it is the fugacity at i? = (tq, -R = ).] 

The function J(t) implements a smearingof the freeze- 
out hypersurface around to; the choice in [B| is 
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(2.11) 



General conditions for J{t) are / dr J{t) — 1, J dr t J{t) 
= To, and At <C to. The last one ensures that one 
can treat fi and H as T-independent, and that freeze- 
out times T < play no physical role. 

For single-hadron spectra J(t) can be immediately in- 
tegrated over, reducing the space-time integral over the 
emis sion function to the Cooper- Frye form (2.1) with 
(2.8). A non-zero duration of particle emission (At > 0) 
has, however, an effect on cluster formation and on 
other two-particle correlations. In two-pion corre- 
lation data from Pb-f Pb collisions at the SPS were fitted 
with At = 1.5 fm/c, although with considerable uncer- 
tainty [Q. This value appears to be small enough to be 
able to neglect the T-dependence of the parameters in fi 
and H; since estimates in have shown that the effect 
of At > on cluster formation should then be small, we 
will continue to use the simpler Cooper-Frye formalism 
(^.3D also for cluster spectra. 



C. Cluster spectra from the model source 



Inserting the expressions from Sec. [I A into Eq. ( |2.3| ) 
one is led to the following integral: 
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The integration over (f) = (j) — ^ can be done analytically, 
yielding a modified Bessel function, but this is not very 
helpful for further analytical progress. More useful is the 
following observation for fireball temperatures 

below the deconfinement temperature of about 150 MeV, 
the ratio Mt/T > 6 A in the exponent is large, and the 
integral will receive contributions only from narrow inter- 
vals T] and p. We may thus expand the hyperbolic cosine 
and sine terms in the exponent, keeping only the leading 
terms: 



ch (ry - r) ch?7t w 1 + i(77 - Y)^ + W 



cos(/)sh?7t ~ rif 



Ap 



(2.13a) 
(2.13b) 



This decouples the longitudinal and transverse flow, and 
the r/-integration can now be done anal ytica lly. Note 
that the "saddle point approximation" ( 2.13| ) becomes 
unreliable for large Pt and/or ijj ||4^ . The remaining 
integrations are easy and give 
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and 
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(2.14b) 



(2.15) 



(2.16) 



In going from ( ^.14a]) t o (|2.14b| ) we replaced the terms 
in the second line on (2.14a) by 1 and assumed non- 
rclati vistic transverse cluster velocities, vt = Pt/Mt ^ 1. 
(Eq. (2.16) should therefore nothe applied to pions |^^!). 
Ves{A^Mt) is the effective volume contributing to the 
emission of clusters with mas s number A and transverse 
mass Mt [0 (see Sect. |l^). T* is the inverse slope 
("effective temperature") of the transverse momentum 
spectrum. 



D. Slope of the Mt-spectrum 



The relation ( 2.16D (with M/A replaced by the hadron 
mass m) was suggested by Nu Xu et al. as a basis 
for a systematic separation of collective flow (77/) from 
thermal motion (T) using single-particle spectra. It has 
recently become very popular, especially since in heavy 
ion collisions at the AGS and SPS the measured m*- 
spectra of most hadronic species seem to follow it quite 
nicely by showing, at sufficiently low pt, inverse slope pa- 
rameters which rise linearly with the hadron rest mass 
g. This has been interpreted as strong evidence for the 
existence of transverse collective flow [ p9t^0|j5l|] . 

When applied to nuclear clusters, however, Eq. (2.1i) 
predicts exactly the same slope for all mass numbers A 
since M ^ A. This is in contradiction with experiments 
both at the AGS and the SPS which show considerably 
flatter mt-slopes for deuterons than for protons [^-p4|. 
In order to understand the possible orig ins o f this dis- 
crepancy one must remember that Eq. ( 2. IE ) rests not 
only on non-relativistic kinematics, but also on the Gaus- 
sian transverse density and linear transverse flow rapidity 
profiles. These were selected in [^,||J|,^,^ for techni- 
cal convenience, and we have kept them here for ease 
of comparability. It was pointed out, however, by Pol- 
leri et al. Q that the Gaussian density profile gives too 
much weight to the center of the fireball where the t rans- 
verse flow is weak. Since according to Eqs. (2.8a) and 
( p.3[ ) the profile function H{p) enters the expression for 
nuclear cluster spectra with the A^^ power, the region of 
weak transverse flow receives the more weight the heavier 
the cluster; for Gaussian profiles this divides the effective 
strength of the transverse flow exactly by a factor A js^ j 
and thus causes the A-independence of the cluster slopes. 

As we will show in Sect. VII B this cannot be com- 
pensated by decreasing the power a in the transverse 
flow profile (2.7) even though this does lead to a more 
rapid increase of the transverse flow rapidity at small 
values of p. Phenomenological consistency with the ob- 
served v4-dependence of the cluster Aft-slopes can only be 
achieved by selecting density proflles H[p) which, when 
taken to the power, don't lose weight in the relevant 
region of large transverse flow. A transverse box profile 
(or a smooth version of it) sat isfies t his requirement and 
is shown to work well in Sect. VII B 



The Gaussian transverse density profile is problematic 
also for a different reason: in combination with the trans- 
verse flow profile ( ^ ) it leads to acausal behaviour for 
heavy particles with large transverse velocities vt- Due 
to the Bolt zma nn factor exp(— P • u{R)/T) in the in- 
tegrand of (2.1) or ( |2.3| ), which gives strong weight to 
space-time points R with u{R) = P/M (the more so 
the larger M/T), such particles are emitted mostly from 
fluid cells at large transverse radii p, i.e. in the tail of 
the Gaussian density distribution. Such matter should 
not exist, however: a cold Pb nucleus with an rms ra- 
dius pi-nis = 5.5 fm |^6|, corresponding to a hard sphere 
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radius po = -y/S/S /Orms, can expand in time tq at most 
to a maximum transverse radius of pmax = Po +to (even 
less, if the transverse expansion velocity is < c). Thus, 
no matter should exist at radii p > Pmax- In our case 
nuclcons with transverse velocities v-t > 0.6 c tend to be 
emitted from causally forbidden regions. (For the lighter 
pions the problem is much less severe due to the larger 
thermal smearing.) 

For the Gaussian density profile ( 2.81| ) we will there- 
fore restrict our attention to nucleons and clusters with 
Vt < 0.6 c, i.e. 1.0 < Mt/M < 1.25. For a proper descrip- 
tion of clusters with larger velocities the Gaussian trans- 
verse density profile must be modified, either by cutting 
it off by hand or by replacing it with a causally consis- 
tent box profile. Unfo rtunately, this forfeits the simple 
analytical expressions (2.14- 2TT6| ) and the direct compa- 
rability with the published results from HBT analyses of 
two-pion correlations [pp|,^ . 

As for the longitudinal rapidity spectrum of clusters, 
we expect from Eq. (2.14) for its width a decrease with 
compared to nuclcons. 



E. Effective source volume Voff and relation to HBT 



The Boltzmann factor in (2.8a) couples the particle 
momentum to the flow vector u{R). This causes a cor- 
relation between the velocity and the spatial coordinates 
of the particle, with a "coupling constant" M/T which 
increases with the particle mass. Particles inside the fire- 
ball are thus sorted with respect to their velocities, and 
particles of given momentum are localized in regions of 
the fireball where the fiow velocity is close to the particle 
velocity. 

Thus only a fraction of the total fireball volume Vcov 
is able to emit particles with given momentum. It is 
this "homogeneity volume" T4iom(''Tit) which is accessi- 
ble through HBT measurements 1^. The HBT radii 
7?.|| (m*) and TZ±{mt) which can be extracted from the 2- 
particle correlation function in the YKP parametrization 
p7| describe the corresponding longitudinal and trans- 
verse lengths of homogeneity in the source. They can be 
evaluated for the model (2.8) as space-time variances of 
the source using the general expressions given in [s^ . If 
these variances are evaluated in the saddle-point approx- 
imation ( |2.13| ) one finds 



^11 = 



Ap 



1 + 

Tq A?7 
Vl + ^(Ar7)2- 



(2.17a) 
(2.17b) 



Thes e are just the factors occurring in the effective vol- 
ume (|2.15|) for A=l: 



Kff(l,mO = (2^)3/2 7^||(m07^i(mt) 
= (2^)3/2 -l4o,„(mt). 



Thus the effective volume T4ff in the cluster spectrum 
( 2.14 ) is very closely related to the homogeneity volume 
extracted from HBT measurements with pairs of identical 
hadrons: 



K=ff(A,Mt) 



A3/2 



3/2 



(2.19) 



For deuterons this implies an effective volume which is 
about 1/3 that of the nucleons. 



Since at AGS and SPS energies A77 > 1 (see Sect. |HF| ), 
the longitudinal flow term ~ {mt/T){A.rjY in TZ\\ dom- 
inates over the geometric term ~ 1. In the transverse 
direction the flow term ^ {mt/T)ri'j is much smaller, 
and the two terms compete with each other, depend- 
ing on the value oi rjf. Higher temperature increases, 
larger transverse flow decreases the lengths of homo- 
geneity. For flxed T, 77/ the homogeneity lengths de- 
crease with increasing transverse mass mt- According 
to Eqs. (2.18,2.19) this implies for fixed particle mo- 
menta a decrease of Vcff (mt) with the particle rest mass 
roughly like M^^^^. Due to the non- negligible geomet- 
ric contribution in TZ± the decrease is actually somewhat 
weaker; for reali stic pa rameters (see below) the combina- 
tion MtVcff in ( ^.14b| ) turns out to be practically inde- 
pendent of Mt: MtVcs{A, Mt) « MVcs{A, M). 

Eqs. (2.17- 2T9| ) indicate perfect mt-scaling of the HBT 
radii and the e ffective volume. If true, the values of VeS 
to be used in (2.14b) for nuclear clusters could be di- 
rectly extracted from two-pion HBT measurements at 
very high pt such that the transverse mass Mt would 
be the same. Unfortunately, this is an artifact of the 
saddle-point approximation (2.13) |4^; a numerical eval- 
uation shows that mt-scaling of the HBT radii is broken 
by transverse flow, albeit weakly [p7l3l. Appropriate care 



(2.18) 



must thus be taken before using (2.18) to compare the 
effective volume Vcff in cluster formation with the homo- 
geneity volume T4ioin extracted from HBT measurements 
with pions or kaons. 



F. Parameters for central Pb+Pb collisions 

We close this section by specifying the model param- 
eters to be used later in the calculation of cluster yields 
and spectra. In Pb-|-Pb collisions at 158 GeV per nu- 
cleon, mid-rapidity is at j/cm = 2.91. For very central col- 
lisions (4-5% of the total inelastic cross section crJot', cor- 
responding to impact parameters 6 < 3.5 fm), the analy- 
sis of pion spectra and two-pion correlations from NA49 
JlO| , ^ 59 1 led to the following estimates for the model 
parameters §,||l: A77 « 1.3, Ap « 7fm, tq « 9fm/c, 
Ar sa 1.5fm/c, rjj w 0.35 and T « 130 MeV. To esti- 
mate the effects of possible errors in T and ryy, we will 
additionally test the {T,r]f) combinations of (100 MeV, 
0.43) and (168 MeV, 0.28), which also describe the slope 
of the single particle wt-spectrum of negatively charged 
hadrons (mainly pions), but not the behaviour of the 
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transverse HBT radius TZ\{mt). The effects of uncer- 
tainties in the determination of tq and Ap, which may 
also be important, have not been studied. From the 
above paramete r va lues we calculate a total covariant 
fireball volume ( |2.9| ) of Vcav ~ 9000 fm'^, homogeneity 
lengths of 7^|| (1 GeV) « (3.2 ± 0.4) fm and 71^(1 GeV) « 
(5.1 ± 0.8) fm, and an effective volume for nucleons of 
V;ff(l,TO) w (0.15 ±0.06) X Kov 



III. THE QUANTUM MECHANICS OF 
COALESCENCE 

A. Reference frames 

So far we have treated nuclear clusters as point-like 
particles without internal structure. This is not realistic: 
the deuteron, for example, has already an rms radius of 
almost 2 fm which is not much smaller than the ho- 
mogeneity radii given above. A proper inclusion of finite 
size effects and the internal cluster wave function in the 
description of the coalescence process is therefore manda- 
tory. This requires a quantum mechanical treatment. In 
this section we will concentrate on the two-body problem, 
i.e. the coalescence of (anti)deuterons. 

Due to its small binding energy, only nucleons with 
small relative momentum can form a deuteron. In the 
deuteron rest frame (d-frame) the relative motion of the 
two nucleons and their coalescence into a bound state are 
thus described by non-relativistic quantum mechanics. A 
relativistic formulation is needed, on the other hand, to 
describe the motion of the deuterons (i.e. of the center 
of mass of the nucleon pair) in the fireball rest frame 
(f- frame). 

In the f-frame, we denote the deuteron's momentum 
and cm. space-time coordinates by Pd = (-EdiJ'd) and 
Rd = (i?j,i2d), and the momenta and space-time coor- 
dinates of its two nucleons by P± and R± = (i?° ,il±). 
In the d-frame, the four-momenta q± of the two nucleons 
are given in terms of their relative momentum q by 



q^,±q 



(3.1) 



Their space-time coordinates r± — (td,r±) in the d- 
frame can be expressed in terms of the d-frame cm. 
coordinates = itd,fd) and their relative distance 
r'^ = (0, r) as r± = rd ± ■^r. Note that in the d-frame 
r has a vanishing zero component, consistent with the 
equal-time nature of non-relativistic quantum mechanics. 
The deuteron's four-velocity in the f-frame we denote by 
6^ = (6°, b) = P^/nid- The matrix L for a boost from 
the f-frame into the d-frame is then given by 



( bo 

-by 





by 


bz 


^ l+fco 




b^b^ 


l+b. 


1+bo 




1+bo 


bijbj 


1+fco 






byb^ 


1 

1+b 


1+fco 


1+bo 



(3.2) 



The matrix L^^ for the inverse boost follows by substi- 
tuting bi ^ —bi for i — x,y, z. 

We consider the f-frame coordinates P± and R± as 
functions of the f-frame cm. coordinates Pd and Rd and 
the d-frame relative coordinates q and r: 



R'^ = {L-r''r^..=K±h{L-'r'' 



P^ = {L-r q±. 



(3.3a) 
(3.3b) 



For the internal (relative) wave function Lpd of the 
deuteron we will consider both the Hulthen form with 
the parameters a=0.23 fm~^, (3—\.Ql fm^^, and the wave 
function of a spherical harmonic oscillator with the size 
parameter (i=3.2fm (r = |r|): 



¥'d(r-) 



a/3(a-f/3) e 



-/3r 



y 27r(a-/3)2 
^d{r)^{iTd?y"' 



exp 



r 

2d? 



(3.4a) 
(3.4b) 



While the first choice gives a better description of the 
deuteron ground state, the second one is technically ad- 
vantageous since it allows for a largely analytic evaluation 
of the coalescence factor below. Both wave functions re- 
produce the measured rms radius of 1.96 fm Please 
note that the variable r describes the diameter and not 
the radius of the relative wave functions such that the 
rms radius is given by r^^^^ — J (Pr (0 |(^d(''")P- 



B. The density matrix formalism 

According to the rules of statistical quantum mechan- 
ics [ pH the number of created deuterons with momentum 
Pd is given by projecting the deuteron density matrix 
onto the two-nucleon density matrix in the fireball at the 
freeze-out time t{: 

^ ^ y"d^xid^a;2 d^x[(fx2(t>*dixi,X2)(l)dix[,X2) 

x{ib\x'2,ti)^^{x[,ti)ip{xi,ti)ib{x2,t{)) . (3.5) 

The correct spin and isospin factors will be added later 
after introducing Wigner functions. A rigorous ansatz 
with projection operators that select the correct spin and 
isospin state for the deuteron from the two-particle den- 
sity matrix can be found in p9| . 

In the detector the deuteron is observed as a free mo- 
mentum eigenstate. Its wave function is therefore given 



7 



by a plane wave for the cm. motfon multiplied by the in- 
ternal wave function (^d : 4'dixi,X2) = (27r)~^/^ exp(iPd- 
{xi+X2)/2) (pd{xi—X2)- The two-nucleon density matrix 
in the fireball is not known and has to be approximated. 
We assume that at freeze-out the nucleons are uncorre- 
lated: 



{lb\x'2,t),P^x[,t)i;{xi,t)i;{x2,t)) 

= {^\x'2,t)i:{x2,t)) {i:\x[,t)^{xut)) 



(3.6) 



The one-particle density matrix can be expressed through 
the 1-body Wigner function as ||6l|,|62| 



{y^Hx',t)i;ix,t)) 



{2n 



cxp (ip • (a;— a;')). (3.7) 



To evalute the integral ( |3.5| ), we introduce new coordi- 
nates r+ = ^{xi + x[), r_ = \{x2 + a^j), i ~ Xi - x\ 
— {x2 — x'r^ and p = \{xi — x'l + X2 — x'2). The first 
two are the classical coordinates of the coalescing nucle- 
ons and will be further reexpressed by and r. For 



the momentum vectors pi and P2 introduced via (3.7) 
we define the relative momentum q = ^{Pi ^ P2) and 
the total momentum P^i = Pi + P2', the identification 
with the deuteron's momentum is due a corresponding S- 
function for 3-momentum conservation which arises from 
the (Pp integration. The d'^^-integration leads to a term 
which we recognize as the Wigner transform of the inter- 
nal deuteron wave function, 

V{r, q)^J d^^ « (r + |) (r" - |) , (3.8) 



which has the following normalization property: 

f -fl-V{r,q)= [ d''rip*,{r)ipd{r) = l. (3.9) 



For the deuteron spectrum we thus obtain 



dNd 3 



d^q d^T 



V{r,q) 



x/p'^(g+,r+)/,r(<Z-,r_), (3.10) 

where q± , r± are functions of the integration variables 
according to Eqs. (3.3). 



C. The problem of energy conservation 



Since in Eq. ( |3.5| ) the two-nucleon density matrix is 
projected on the density matrix of an energy eigenstate 
(the free deuteron), energy is clearly conserved. On the 
other hand, if we were to replace the Wigner functions in 
( |3.10 ) by classical distribution functions of on-shell nu- 
cleons (as done e.g. in j2l]) we would violate energy con- 
servation: due to the deuteron binding energy two free 



nucleons cannot coalesce into a deuteron without help of 
a third body. In o ther words, one or both of the Wigner 
functions in ( 3.10D are probed off-shell, and the nucleons 
in this expression cannot be considered as free particles. 
(For free distinguishable particles in a thermalized sys- 
tem the Wigner function is identical to the Boltzmann 
distribution ||6l[|.) 

A formalism which makes the conservation of energy 
during coalescence more explicit by including the scat- 
tering with a third body and which thereby gives more 
control over the replacement of the Wigner functions by 
classical distributions was developed by Danielewicz and 
Schuck p^. Instead of (|3.5|) they start from p3| 



— — = lim 



1 1 



TT' 2! 



2T 



:-T / dt 



2T' 



dt' / d^xi d X2 



X d^x[ d^x'2 e''^^^'-*'^(l)*i{xi,X2) Mx'i,x'2) 
y.{i!\x'2,t')i!\x\,t')i!{xi,t)i!{x2,t)) . (3.11) 

They then consider the equations of motion for the 
correlator l^ij]'^ {x'2,t') {x'i,t') ^{xi,t) 'ili{x2,t)Y These 
include interactions with third particles which can put 
one of the two nucleons slightly off-shell. This nucleon 
can then form a deuteron with an on-shell nucleon with- 
out violating energy conservation. It is convenient at this 
point to introduce a four-vector notation, with momen- 
tum four-vectors which can be either off-shell (indicated 
by an asterisk) or on-shell (no asterisk): pi = {Ei,pi), 
p* — {uj,pi), and ri = {ti,ri). For the evaluation of 
(|3.1lD we refer to Ref. [||; the result is (see Eqs. (16,19) 



dNi _ -3i 
~ (27r)3 
x(27r)4(5''(Pc 



dVd d^r 



d'^Pi d^p2 



(27r)4 (27r)3 



-p\-p2)V{r, 



P1-P2 ' 
2 / 



(3.12) 



E<(pt, r+) f^{p2, r_) + E<(pJ, r+) f^{p2,r^)) , 



where the terms ^^f^ and give equal contribu- 

tions. The rate E^(p*,a;) for producing a nucleon with 
off-shell four-momentum p* at point x is given by the 
scattering of the nucleon with a third particle: 



-iT,^{p*,x) 



E 



d^q d^p' d^q' 



(27r)3 (27r)3 (27r)3 



X (2nf5\p*+q-p'-q') |AfNj^N,f 
X f^{p',x)ff{q',x){l±ff{q,x)). 



(3.13) 



Ref. [gSj considers a gas of nucleons only, and the refor e 
there is no summation over particle species j as in (3.13). 
Our fireball, however, mainly consists of pions, and all 
particle species with large cross sections with nucleons 
must be accounted for in the production rate 



With the 4-dimensional ^- functions in (3.12,3.13) now 
explicitly accounting also for energy conservation, we 
proceed to approximate the Wigner distributions fj^ in 



8 



these e qua tions by the classical thermal distribution func- 
tions ( p.8[ ). Using the momentum-conserving (5-function 
5'^{p*+q—p'—q') we then have the identity 



■■■P'/T 



1± 



1 



-u-p'/T 



1± 



1 



and hence 

-iJ:<{p*,x) ^ fNip*,x)Y^ J -^^fjii^x) 
(Pp' (Pq' 



(3.14) 



x|AfNj^Njf (l±/j(g',.x)) 



(3.15) 



The term in square brackets is the scattering cross section 
times the relative velocity between (off-shell) nucleons of 
momentum p and particles j with momentum q, averaged 
over the thermal distributions: [. . .] = {cFnjV^j{p,q)). 
The remaining integral over dPq and summation over j 
then gives the total scattering rate or inverse scattering 
time of a nuclcon with momentum p at point x: 



(3.16) 



The production rate S^(p*, R) is therefore just given by 
the (off-shell!) nucleon distribution function /N(p*,i?) 
divided by the scattering time of nucleons. 

A deuteron has twice the scattering rate of its con- 
stituent nucleons. Since a scattering event is likely to 
break up the deuteron, and since we have by assumption 
free streaming particles after freeze-out at t{, the time 



integration dt^ in (3.12) should start at t{ — and 



end at t{. Further assuming that we may treat the distri- 
bution functions as constant over this time interval, the 



inverse scattering time in (3.16) cancels against the time 
integration, and wc arc left with 



(27r)3 



xfpiq 



d^r d^q 
(27r)3 

r+) fn{q*-;'i 



V{r,q) 



(3.17) 



The difference between this expression and ( 3.10| ) is that 
the Wigner distributions have been replaced by classical 
distribution functions (2.8) one of which is off-shell (gl) 
so that a deuteron may be formed without violating en- 
ergy conservation. The energy component of g* follows 
from the condition + q* = (md,0) implied by the 
(5- function in ( ^.12 ): 



q1 = {^\/m^ + , 
q_^ = l^md - \/m? + , -qj 



(3.18) 



Eq. (3.17) is also a good basis for studying the quan- 
titative effects of energy conservation on the coalescence 
process, for example by replacing the off-shell momentum 
ql by its on-shcU limit. 



D. Coalescence in the highly relativistic fireball 

We now return to the relativ istic moti on of the par- 
ticles in the fireball frame. In ( [3.10 ) or ( |3.17|) we need 



the arguments of /(p, x) in the f- frame, so we substi- 
tute r± and q± by R± and P± using ( |3.3| ). After per- 
forming the d^r and S'q integrations in the d-frame, we 
reexpress the S't^ integration over the deuteron's cm. 
coordinate at time tf in the d-frame in covariant form: 
E^d^rA = Pd ■ rf^cr(i?d)- The integral d^a{Rd) extends 
over the freeze-out surface Sf which is characterized by a 
relation i?d f(-Rd) between the deuteron freeze-out time 
and point in the f-frame. In this way we obtain 



diVd _ 
d^Pd (2^)3 

X /p (^d, f ) /n (i?d, f) Cd(i?d, Pd) (3.19) 



with 



d3gd3 

72^ 



■Vir,q) 



/(i?+,P+)/(i?„,P_) 



(3.20) 



Up to the quantum mechanical correction factor Cd this 
is identical with the classical formula (2.3). Cd(i?d,Pd) 
is an integral over the internal phase-space coordinates 
and provides a measure for the homogeneity of the nu- 
cleon phase-space around the deuteron cm. coordinates 
(i?d,Pd/2). For a homogeneous nucleon phase-space 
(static and very l arge systems) the second factor under 
the integral ( 3.20| ) (which we will call "homogeneity fac- 
tor") approaches unity and ( [3.2C| ) reduces to the normal- 
ization integral (|3.9[), yielding Cd = 1. 

As we will see, the g-dependence of t he di stribution 
functions f{R+,P+) and f{R-,P-) in ( ^ ) is much 
weaker than that of the deuteron Wigner density I?(r, q) 
which is peaked near q = 0. The distribution functions 
can thus be pulled outside the g-integral, yielding 



Cd(-Rd, Pd) 



/(P+,f)/(P_,f 



P {Rd 



3i 

2 



\Mr)\' 



(3.21) 



This expression is now very similar to the Hag edorn 
model for cluster production in pp collisions 137]] . Ac- 
cording to Hagedorn, the probability of finding a cluster 
with certain quantum numbers (mass M, spin, etc.) is 
equal to the probability to find a really elementary parti- 
cle of the same quantum numbers, times the probability 
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C that the cluster is contained inside the reaction volume 

n-. 



\ip\'dV. 



(3.22) 



For pp collisions Hagedorn assumed that the reac- 
tion volume was given by the pion Compton wavelength, 
= 47r/(3m^), resulting in C ~ 0.17 for the deuteron 
and C « 0.28 for '^He or ■^H. This is consistent with the 
homogeneity radii extracted from two-pion HBT correla- 
tions in pp collisions which are also of the order of the 
pion Compton wavelength KJ . 



IV. THE QUANTUM MECHANICAL 
CORRECTION FACTOR 

In this section we evaluate the coalescence factor Cd 
for the source model discussed in Sec. || and relate it to 
the homogeneity radii TZi which can be extracted from 
HBT measurements. 



A. Integration over relative momenta 

The explicit dependence of the coordinates R± in 
( |3.20 ) on r is given by (3.3): 



From _R-|- one must calculate 



b ■ r 

1 + W 



1 Rl + R±, 
= o 1^ ■ 



2 Rl- R±, 



R±x + 



±y : 



P± 



(4.1) 

(4.2a) 

(4.2b) 
(4.2c) 



which enter the density profiles H{R±) and th e flo w vec- 
tors u{R±). Since the Lorentz transformation (3.3) mixes 
the components of r, the fu nctio ns ?7±(r'), p±{r), T±{r), 
and thus the integrand of (3.2C) are in general compli- 



cated functions of the integration variables r^, Vy und 
rz- 

The argument u{R±) ■ P± of the Boltzmann factor can 
be evaluated in any frame. It turns out to be most con- 
venient to transform u{R±) and u{Rd) from the f-frame 
to u± and in the d-frame: 



= U-'u.iR^), < = L^^u.iRi). 



(4.3) 



Under the assumption that over the effective integration 
domain the flow and chemical potential can be approxi- 
mated as constants we find 



/(i?+,P+)/(i?_,P_ 



P (i?d. 



exp 



X exp 



T 



2(Ap)2 



+ ri-- H 
2(A77)2 



(4.4) 



Since in the d-frame the relative momentum q is small, 
the energy components q^. can be expanded non-relativi- 
stically: 



q+-u+ + q^-u^ = m(u°, + u°_) + (w° ± w° ) 

2m 

-q-{u+-u_) + ^{u^_Tu'L), (4.5) 

where B — — 2m. The two signs refer to two different 
treatments of energy conservation: the lower sign applies 
if one of the tw o nuc leons (here: g_) is off-shell as pre- 
scribed by Eq. ( 3.18 ), the upper sign refers to the case 
where we simply neglect energy c ons ervation and take 
both nucleons on-shell as in Eq. (^]^) (in this case the 
term ~ B vanishes) . With harmonic oscillator wave func- 
tions the Wigner density of the deuteron is a Gaussian, 



T^i'i'iO) = 8 exp(— r /d — q c? ), and the integration over 



(Pq in (3.2C) is straightforward: 



Cd(i?d, Pd) 



dh 



P\ + P- - 2p 



2(Ap)2 



X exp 



X exp 



r exp 

2(Ary)2 

,.2 

III . 

Bcs_ {u+ - u 
T 



xexp(-^-^«+«" -2.°; 



iT^dl 



ff 



(4.6) 



Here the two parameters d^^ = d"^ + {u^ ± vP_)/{2mT) 
and Bcs — B{vP_ =F )/2 depend again on how we deal 
with energy conservation. 



B. Integration over relative coordinates 

When trying to perform the integration over d^r one 
encounters subtle causality pro blem s which require some 
discussion. The integration in ( 3.19] ) runs over the freeze- 
out hypersurface Td = tq. The hypersurface spanned by 
the d^r integration at constant time td in the d-frame is 
given in the f-frame by 



{R-RA)-blb\ 



(4.7) 



which follows by substituting b R± into i?^ in Eq. ([l.l|). 
This hypersurface is inclined relative to the freeze-out 
hypersurface, and it also cuts the T = \Rz\ half-planes 
which are tangent to the light cone originating from the 



10 



collision point of t he tw o nuclei. Therefore, when inte- 



grating over d r in ( 3.2C ) looking for nucleons which may 



form a deuteron, the coordinates R± leave the freeze-out 
hypersurface and even the light cone which contains all 
produced particles. 

It turns out that this problem is much less severe than 
it first appears. First, contributions near the longitudi- 
nal light cone are suppressed by the fact that on the light 
cone the longitudinal proper time r = (i?'')^ — va- 
nishes, leading to a diverging flow velocity u whose tem- 
poral and longitudinal components are given by R^ /t and 
Rz/t, respectively, and, correspondingly, to a vanishing 
Boltzmann factor exp(—p-u{x)/T). This constrains the 
effective integration domain in the longitudinal direction 
to a region well inside the light cone. 

In the transverse directions we must constrain the in- 
tegration domain by hand to th e inside of the light cone, 
since (as discussed in Sec. [ID) our source parametriza- 
tion does not automatically ensure consistency with 
causality of the produced particles' positions and veloci- 
ties. Fortunately, this is not a serious interference: both 
the deuteron Wigner density and the "homogeneity fac- 
tor" in (3.2C) are peaked at small values of r. The lat- 
ter is the more peaked the more strongly the system ex- 
pands, i.e. the larger the flow velocity gradients are. As 
we will see, the integration domain over r is, to the ex- 
tent that it is not already restricted by the finite range 
of I?(r, q), aga in gi ven by the homogeneity radii TZi dis- 



cussed in Sec. HE. Fortunately, these are smallest for 



deuterons with large transverse masses where the causal- 
ity problems are expected to be most serious. 

In view of their minor numerical effects, we refrain 
from giving a detailed account of the technical implemen- 
tation of these restrictions into the numerical evaluation 
of (4.6), referring instead to Ref. [Q. In practice, the 
following analytical estimates turn out to be sufficiently 
accurate even on a quantitative level. 



C. Analytic approximation of the correction factor 

Since the measured deuteron momentum spectra do 
not contain information on the point of deuteron forma- 
tion, the relevant quantity is the average correction factor 



(Cd)( 



p . ^ J Pd-rf'g(fld) P {Rd, ^) Cd(i?d, Pc 

/Pd-dM^d)/M^d,f ) 



(4.8) 



We will first calculate Cd(-Rd, ^d) in analytic approxima- 
tion for a special combination of i?d and Pd (given by 
(4.9) below) and then argue that the result, called Cj, is 
actually a very good approximation of (Cd) (Pd)- Numer- 
ical studies |4^ confirm the validity of these arguments. 

We first concentrate on deuterons with zero transverse 
momentum which are at rest in the fluid cell where they 



are created, i.e. whose four-velocity 6'' agrees with the 
flow four- velocity at the production point, b — u(i?d): 



6^ = -(i?J], 0, 0, i?d.) 



(4.9) 



In the d-frame the fireball near i?d can then be described 
non-relativistically as long as the longitudinal flow veloc- 
ity gradients are sufficiently small. As long as d < Td the 
non-relativistic approximation is very good in the rele- 
vant region d where T) is non-zero. Since bx^by~0 



the Lorentz transformations (3.3,4.3) do not mix longi- 
tudinal and transverse directions, and the integrand for 
the coalescence factor has the same axial symmetry as 
the fireball. Then 



7/± W 7/d ± — , 

r'^ -I- 
2 2 , ^ V 

p± = Pd + — : — 

T± W Td 



± {rxRdx + ryRdy) , 



8Td 



(4.10a) 
(4.10b) 
(4.10c) 



after a non-relativistic expansion of ( |4.2| ). Further, the 
flow u± follows from u{R±) by a simple shift in the lon- 
gitudinal rapidity: 



iv-Yd? 



9 9 

VfP 



VfRx VfRy 

'2(Ap)2' Ap ' Ap 



where we have already used the saddle-point approxima- 
tion ( ^.13|). For the given values of d and mT the value of 
doff in ( |4.6| ) depends only weakly on Since Uj.(r)>l 
for b—u{Rd) and small |r|, we use d'^^^d'^ +1 / mT and 
-Scff=0, or d^ffKid'^ and BcS~B, respectively, depend- 
ing on whether or not energy conserva tion is taken into 
account. With these approximations ( |4.6| ) turns into a 
product of Gaussian integrals in Tx, Vy and r^, with the 
result 



-L / " \ / ^cff 

Cd = — I ^— I exp ( 



71711 V4ff 



7± - W 1 



711 = 




2n±{m)J \2TAp 



1 



27^||(m)y \2TTd 



(4.11a) 
(4.11b) 
(4.11c) 



For the source parameters given in Sec. |IIF| the last terms 
under the square root in 7|| are negligible . T hey orig- 
inate from the coupling term g - (m ^— tt_) in (4^) and the 
resulting term it_p in (4.6). They would th us be 
absent if we had started from the approximation (3.21) 
instead of ( 3.20| ). The sm allnc ss of these terms is a good 
check of the accuracy of ( [3.2l[). 

The last two factors in (4.11a) deviate from unity by 
less than 2% for temperatures T between 100 MeV and 
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170 MeV if energy conservation is properly accounted for; 
if not, the deviations have the opposite sign, but remain 
below 5%. 



Given the high accuracy of ( 3.21 ), we can use it for 
a check of the sensitivity of C° on the choice of the in- 
ternal deuteron wave function. A numerical integration 



of (3.21) with the Hulthen wave function (3.4a) yields 
values for which are somewhat larger than those for 
harmonic oscillator wave functions. For the source pa- 



cillator wave function 



rameters given in Sec. II F we obtain for the harmonic os- 

+0.03 



form 



0.84 



+0.02 
-0.04 



0.811^:^^ and for the Huhhen 



(where the upper and lower limits 
indicate the effects from the estimated uncertainties in 
(T, r//)). The numbers in Table | show that the differ- 
ences are sensitive mainly to the transverse and longitu- 
dinal flow gradients; they remain on the level of a few per- 
cent for weakly expanding sources (leftmost column), be- 
come stronger for more rapidly expanding sources (right- 
most column) , and can reach a factor of 2 or 3 in systems 
with very small interaction volume {pp collisions). 

The origin of the difference is readily understood: 
while both wave functions provide the same rms radius, 

1.5 fni for the 
3 fm for the harmonic oscil- 



- C° in front of that infe- 
ct. A numerical check gave 



the maximum of r |iy9d('')| 

Hulthen form and at ', 

lator. Since the "homogeneity factor" in (3.21) peaks at 
small values of r, especially for strongly expanding sys- 
tems with small homogeneity radii, the integral is larger 
for the more realistic Hulthen wave function than for the 
harmonic oscillator one. 

Numerical calculations show that Cd{Rd,Pd) varies 
much less as a function of than /^(i?d, -Pd/2). On 
the other hand, the particular point with w(_Rd)=& 
at which was evaluated corresponds to the maximum 
of Cd(i?d,^'d), to the maximum of the Boltzmann part 
of /(i?d, ft/S), and thus approximately to the maxi- 
mum of integrand in the numerator of (4.8). We can 
therefore pull Cd{Rd,Pd) 
gral and thus have (Cd) - 
(Cd) (Pd=0) = 0.79 instead of CJJ = 0.81. 

For a boost-invariant source, where every comoving ob- 
server has identical surroundings, we expect Cd(i?d,^'d) 
to depend only on the difference between the local 
flow velocity and the dcuteron's velocity, Cd{Rd:Pd) = 
Cd{u{Rd) — b), and (Cd) to be independent of Pd- In 
our fireball longitudinal boost-invariance is broken by the 
density profile H{R). However, for the Gaussian profile 
used here (Cd) still turns out to be independent of the 
deuteron's longitudinal rapidity. In the transverse direc- 
tion we have no boost-invariance at all. As a consequence 
we find a slight decrease of (Cd) with increasing trans- 
verse velocity of the deuteron. Deuterons with non-zero 
transverse velocity see the fireball Lorentz-contracted in 
their direction of motion; this decreases the correspond- 
ing length of homogeneity and thus (Cd). In the region 
nit/m < 1.25, i.e. Vt < 0.6 c, to which we restrict our dis- 
cussion in the case of Gaussian transv erse d ensity profiles 
(see the discussion at the end of Sec. IIP), this effect is 



small and (Cd) is approximately constant. 

We can summarize the results of this section in the fol- 
lowing approximate formula for the quantum mechanical 
correction factor in terms of the deuteron size parameter 
d and the longitudinal and transverse lengths of homo- 
geneity for nucleons: 



1 




(Cd 



(4.12) 

This expression does not depend on the longitudinal ra- 
pidity and, for small transverse velocities, only weakly 
on the transverse momentum of the deuteron. Apply- 
ing this expression to p-l-p collisions and inserting corre- 
spondingly for the homogeneity lengths about 1 fm each 
one obtains (Cd) ~ 0. 15 in good numerical agreement 
with Hagedorn's value ( 3.22 ). 



V. LARGER CLUSTERS: ?H AND ^HE 

In this section we give a quick estimate of the quantum 
mechanical correction factor for clusters made of three 
nucleons, i.e. for '^H and ^He. Clearly, for three nucleons 
the internal wave function and its Wigner transform is 
much more complicated than for deuterons. For the three 
nucleons with coordinates /"i, T2 and r^, we introduce the 
cm. coordinates R = (^i+f 2+''"3)/3 and the relative co- 
ordinates p = (ri— r2)/\/2 and A — (ri-f r2— 2r3)/\/6. 
With this choice we have r\+r2+r\ = R'^+p^+X^ and 
d'^rid'^r2d^r3 = 3^ d'^Rd'^pd^X. As before we approxi- 
mate the internal wave function by a spherical harmonic 
oscillator solution M: 



'Piri,r2,r3) (3T:'^b*) ^^"exp 



2&2 



(5.1) 



This wave function is normalized and has the rms radius 
b: 



33/2d3pd3A|(^(ri,r2,r3)P = 1, 



Both clusters are spin-^ fermions, and the binding en- 
ergies and rms radii are approximately —8 MeV and 
1.75 fm. (Note that this rms radius is smaller than for 
deuterons!) '^He is somewhat more loosely bound than 
■^H, but we neglect this difference here. 

To estimate C° for a "^H/^He cluster at rest in the center 
of the fireball, with momentum P = (3m, P—0) and cm. 
coordinates R = (i?",i?=0), we evaluated numerically 
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f{R,,i)f{R,,i)f{R„i) 



(5.3) 



where Ri = (i?", r^) are the space-time coordmates of the 



three nucleons. The results C'^=0.78to°Ql for To=9fm/c 
and Ct=0.67lQ 07 for ro=6fm/c, respectively, are not 
much smaller than the corresponding values for the 
deuteron. Whereas one would generically expect larger 
flow effects for three than for two nucleons, the three 
nucleons in the '^H/^He clusters occupy a smaller and 
therefore more homogeneous region around their center 
of mass. 

As for the deuteron we expect that C° provides a good 
estimate for the average correction factor (Ct). A more 
rigorous calculation must, however, take into account the 
binding energy which is larger and thus more important 
than that of the deuteron. Unlike deuterons, tritons and 
■^He can be formed via an excited state by coalescence of 
three nucleons without requiring additional particles for 
energy-momentum conservation. 



completely in terms of the deuteron size d and the ho- 
mogeneity lengths ("HBT radii") 7^J^, 7^||. 

Eq. ( |6.3| ) implies that, in the model of Sec. Ba 
is almost momentum independent: both mtVcff (1, m,t) 
and (Ca) depend only very weakly on rrit and Y. An 
important precondition for this weak mj-dependence of 
B2 is, of course, the cancella tion of the exp(— Mt/T*)- 
factors; as discussed in Sec. |ll El| , the latter is due to 
the Gaussian form of the tran sverse density profile in 
H{R) which according to ( 2.16 ) causes identical inverse 
slope parameters T* for all clusters. We have mentioned 
before that this is inconsistent with the measurements 



(see Sec. VII B below), and that more box-like transverse 
density p rofil es are phenomenologically preferred. In this 
case Eq. (6^) must be amended as follows: 



Bo = 



3 7r3/2 ^c^^ 



2(mt — m) 



2mt 7^i(mt)7^^|(mt) 



(6.4) 



Eqs. (412,^^4) are the most important theoretical results 
of the present paper. 

Since neutrons are hard to measure, experiments usu- 
ally do not provide Ba, but rather 



VI. EXTRACTING PHYSICS FROM MEASURED 
CLUSTER SPECTRA 

A. The invariant coalescence factor Ba 



Accordi ng to ( [3.19| ), the invariant cluster spectra are 
given by (2.14) multiplied by the quantum mechanical 
correction factor {Ca){P)- 



dNA ^ ^. 2Ja+1 



(27r)3 



{Ca){P) K=ff(A,Mt) 



X exp 



Mt - M 



(6.1) 



For the invariant coalescence factor Ba defined by (1.1) 
we thus find 



2^ 



(27r) 



Vcsil, mt) \mtVes{l,mt) 



A-l 



(6.2) 



The factor A arises from Mt/mf — A/mf^^. (In case 
of a static, non-expanding fireball the homogeneity vol- 
ume VcS in this expression would be replaced by the 



total fireball volume Kov il^ 



Eqs. (2.18, Ogj ), we can write B2 as 
3^3/2 



With Vcff given by 



Bo 



2mt7^i(mt)7^||(m^) 



(6.3) 



Note that the last exponential factor in (6.1), which de- 
pends strongly on Mt a nd Y , has cancelled in the ra- 
tio. With (Cd) given by ( 4.12| ), B2 can thus be expressed 



b: 



dNA 
d^PA 



^dN, 



Ba exp 



"d^Pp 
T 



Z+N 



P^=Pa/A 



(6.5) 



Here possibly different chemical potentials for neutrons 
and protons are important: If nucleons and antinucleons 
have the same temperature, flow and freeze-out density 
distribution, as we have assumed, our model yields iden- 
tical coalescence factors Ba and Ba for clusters made of 
matter and antimatter. For /Xn ^ /^p the corresponding 
values B\ and B*^ will, however, be different, and they 
will also differ from Ba- If the initial neutron excess 
of the cold Pb nuclei were still present at freeze-out, we 
would expect for Pb+Pb coUisions B*^ w 1.5 B^ w 2.3 B^, 
i.e. quite large differences. Of course, it is not likely that 
the large net isospin remains in the neutron channel until 
freeze-out; a considerable fraction is expected to boil off 
with other produced particles. Nevertheless, there may 
be a visible effect of 7^ A*p on P2 ^nd especially on B^ 
for (which contains 2 neutrons), as well as a charac- 
teristic difference in the B3 values for '^H and ■^He. 



B. Cluster fugacities 

From the ratio of particles to antiparticles the fugaci- 
ties can be calculated: 



E 



_ dNA 



exp 



Pa=Pa 



(^) 



(6.6) 
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Using the proton fugacity Ap = exp(/Zp/T) from the ra- 
tio p/p and the deuteron fugacity Ad = exp ((yLtp-|-yLtn)/r) 
from the ratio d/d one can in principle extract the neu- 
tron fugacity An = exp(/in/T). Once the temperature is 
known (see below), the separate chemical potentials of 
neutrons and protons can thus be determined. In prac- 
tice, however, the uncertainties in the p and p spectra 
from insufhciently well known A and A decay contami- 
nations and the large statistical error bars on the d spec- 
trum limit the usefulness of such an analysis. 



C. Freeze-out temperature from cluster ratios 

The thermal freeze-out temperature for pions was de- 
termined from a simultaneous analysis of t heir spectra 
and two-particle correlations, see Sec. 



[IF. The same 



freeze-out parameters appear to also describe quite well 



the proton mt-spectra |66 . An independent determina- 
tion of the proton freeze-out temperature uses the chemi- 
cal composition rather than the shape of the momentum 
distribution at freeze-out. In particular one can try to 
analyse the ratio 



Saa' 



dNA 



Ea' 



dNA' 



(6.7) 



,=0 



of the invariant spectra at zero momentum of two dif- 
ferent types of clusters (including nucleons, ^=1). To 
investigate such a possibility is suggested in particular 
by the data from the NA52 experiment at CERN Q 
and from the E864 experiment at the AGS which 
mea sure cluster yields only very close to Pj = 0. With 
(3.1) this ratio is given by 



Saa' — 



2Ja + 1 



X^A-Z^, X^^-^A' f,iA'-A)^ 

AVMMa) 



2Ja' 

(Ca) 

{Ca') A'V,s{A',Ma') 



(6.8) 



which for known fugacities An,p is easily solved for the 
temperature T. [For antiparticles Z and N must be taken 
negative in (|6.8|).] 

The main error of such a determination of the tem- 
perature does not arise from the correction factors (Ca), 
which do not vary much in our parameter range, but 
from the sometimes substantial uncertainties of the ex- 
perimental value of Saa' and of the fugacities. Further- 
more, the temperature value determined from the yield 
ratio at a certain point in momentum space rather than 
from the ratio of total yields depends on model assump- 
tions about the shape of the spectra. For example, we 
already mentioned that a box profile for the transverse 
density distribution gives better results for the slope sys- 
tematics of the transvers e clu ster spectra than the Gaus- 
sian profile from which (6.1) was derived. (This change 
does not affect the rapidity spectra.) Adjusting the box 
radius such that the same rms radius Ap is reproduced 



(Pbox = 2 Ap) , we find instead of ( |6.8| ) the same expres- 
sion multiplied by a factor ka/ha' where 



KA 



e — 1 ( 1 



T 



(6.9) 



(We used the sam e saddle-point approximation as in the 



derivation of ( 6.1 



see Sec. ||. We also kept the strength 
of the flow r]f/Ap fixed when switching from the Gaus- 
sian to the box profile.) For sufficiently large trans- 
verse flow ^ T/m this ka is just a factor 1 with- 
out any consequences. In the opposite limit, however, 
lim^j^^o Ka = 2 A, i.e. the intercept of the invariant spec- 
trum at Pt = is changed by a factor 2A. Similarly, 
changing the Gaussian longitudinal density distribution 
into a box doesn't matter much for systems with strong 
longitudinal expansion, but yields another factor ^/A for 
static systems. For static syste ms this just compensates 
the factor A~^/^ from Eq. ( 2.19 ), which is intuitively cor- 
rect since for static fireballs with constant density the 
effective volume Vcs must coincide with the total fireball 
volume, independent of A. 

Parametrizing the last factors in Saa' as 



KA AV^s{A,Ma) 
KA' A'V,siA\MA') 



(6.10) 



we thus find find — ^ J; X < 1, with the upper limit repro- 
ducing the behaviour for a static, homogeneous fireball 
and the lower limit corresponding to rapidly expanding 
systems or systems with a Gaussian transverse dens ity 
profile. For the source parameters given in Sec. II F we 



expect to be close to x = However, in order to allow 
for this kind of model-dependence, freeze-out tempera- 
tures extracted from Saa' should be plotted against the 
power X- 



VII. COMPARISON WITH EXPERIMENT 

In this section we discuss the presently available cluster 
data from Pb+Pb collisions at the CERN-SPS, taken by 
the NA52 and NA44 collaborations. The NA52 exper- 
iment does not trigger on collision centrality and mea- 
sures particles close to Pt = 0, with different rigidities 
P/Z of the spectrometer magnets providing data at var- 
ious rapidities. Recently, the minimum bias data were 
re-analyzed to extract the centrality of the collisions. Al- 
though the minimum bias data contain information on 
p, d, "^H, "^He, p, d and •^He their analysis in the 
framework of the present model is difficult: the assumed 
axial symmetry of the source applies only to central col- 
lisions, and also the model parameters given in Sec. II F 
were extracted from very central Pb-|-Pb collisions. For 
this reason only the impact parameter selected data from 
the 4% most central Pb+Pb collisions will be really use- 
ful for us. At the Quark Matter '97 conference prelim- 
inary A, A-corrected particle ratios p/p, d/d, i?2(d/p^) 
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and B2(d/p^) were presented for y = 3.75 as a function 
of centrality jj]. They were recently amended by a cor- 
rection to the proton spectra from A-decays |^^. 

From NA44 we have the Mt spectra of p, d and 
i32(d/p^) as well as an Mfintegrated ratio d/d w 
(0.4 - 1.5) • 10-3 [0 in the rapidity range 1.9 < y < 2.3 
for the 20% most central collisions. 

Finally, we have data for B2 and from Au+Au colli- 
sions at the AGS with similar centrality (4% CTtot), but for 
the considerably lower beam energy of 11.5 GeV/nucleon 
p2[ . Comparison with these data requires a readjustment 
of the source parameters. 

A full comparison between theory and data is possible 
only for the cluster spectra taken in very central (about 
4% (Ttot ) Pb +Pb collisions, since the model parameter set 
of Sec. II F describes this class of collisions only. The 



one data point at this centrality provided by NA52 can 
therefore only be used to test the absolute normalization 
of (6J_) or the value (^|^) of Ba and does not give in- 
formation on the spec tral behavior. On the other hand, 
the Mt-shape of ( pj| ) can be tested qualitatively with 
the NA44 data although, due to the weaker centrality 
cut, the normalization and the slope parameters may be 
somewhat different from what we expect on the basis of 
the parameters from Sec. [IF. 



A. NA52 data 



From the parameters given in section II F and from our 
values (Cd) ~0.8 and (Ct) w (Cshc) ~0.7 for the quantum 
mechanical correction factors, we estimate the invariant 
coalescence factors Ba in 4%crtot Pb-|-Pb coUisions by 
using (6.2): 



B2 
B3 



'+1 

(2;1 



(6;if) • 10"^ GeV^ 



10"^ GeV* 



(7.1a) 
(7.1b) 



Fig. shows the preliminary NA52 data [|| for B2 as 
a function of collision centrality, with the most central 
collisions on the left. The agreement with our estimate 
is very good, both for d and d, but we should note 
that it deteriorates somewhat if the recently reported A- 
corrected preliminary data ||6^ are used for the compar- 
ison. A small neutron excess at nucleon freeze-out may 
explain the sytematic ally slightly lower B2 values for an- 
tideuterons (see Eq. (6.5)). Fig. |^ shows an increase of 
B2 by up to a factor 8 for more peripheral collisions, 
in qualitative agreement with the naive expectation of 
decreasing homogeneity lengths and a smaller effective 
volume in less central collisions. 

From the preliminary p/p and d/d ratios at y = 3.75 
B and a preliminary d/p ratio |^ in the most central 
(4% (Ttot ) impact parameter bi n we tried to determine the 
freeze-out temperature using (6^). If we assume a com- 
mon nucleon fugacity corresponding to 11 /T ~ 1.5, we 
obtain T « (144 ± 5) MeV for Gaussian density profiles 



{X= - 0.5) and T « (124 ± 4) MeV for a weakly expand- 
ing homogeneous fireball (x=lj but still with (Cd) = 0.8). 
For a larger nucleon chemical potential of /x/T=1.9 we 
find T « (135 ± 5) MeV and T « (118 ± 4) MeV, respec- 
tively. The given range of fi/T at thermal nucleon freeze- 
out and the resulting uncertainties for the temperature 
are due to the experimental errors on the particle ratios. 
Please note that the different values for x cause a differ- 
ence of almost 20 MeV in the extracted temperature. For 
a purely thermal model with x = 1 ^nd (Cd) = 1, NA52 
give a freeze-out temperature of T « (115 ± 10) MeV 
H, taking all uncertainties into account. For a trans- 
verse box and longitudinal Gaussian density profile (as 
will be motivated below), we have x ~ 0-5, and using 
again (Cd) = 0.8 we obtain T « (125 ± 10) MeV. Given 
the systematic uncertainties, all these values are consis- 



tent with the model parameters given in Sec. II F 



B. NA44 data 

In Fig. 1^ we show the Mt distributions of protons, 
deuterons, and B2 from our model. These figures grew 
out of an analysis of preliminary data by the NA44 col- 
laboration which were made available to us by M. Murray 
1^ , but have so far not been published (which is the rea- 
son why no data are shown) . A fit of the Mj-dependence 
of the invariant proton and deuteron momentum spectra 
with the function exp(-(Mt - M)/T*) yielded effective 
temperatures T* « 250 MeV and « 350 MeV, respec- 
tively, i.e. a considerably larger value for deuterons than 
for protons. As a consequence, the measured B2 rises 
with Mt by about a factor 2-3 in the transverse mass re- 
gion mt/m =1.00-1.25. Higher values of T* for deuterons 
than for protons have been repeatedly observed p^-|5^] 
and are in clear contradiction with the prediction of iden- 
tical slope parameters (2.16) from a Gaussian source 
model. 

Recently, several possibilities were suggested to ac- 
count for this discrepancy | |55| : a transverse flow profile 
(^.7|) with a = 0.5, a box profile for the transverse den- 
sity distribution, and surface emission of clusters. We 
have analyzed the first two possibilities. 

Starting from (2.12) with a linear transverse flow pro- 
file a = 1 and the source parameters of Sec. |IIF| , we 
first selected different combinations of (T, rjf, ^ which 
all gave good fits to the proton spectrum (Fig. ^). The 
lines for the deuteron and B2 spectra in Figs. ^,c follow 
then self-consistently from the parameters used in the 
proton spectrum without further adjustment. Line no. 6 
in Figs. ||b,c corresponds to the parameter triplet (140 
MeV, 0.34, 1.75); it clearly fails to reproduce the deuteron 
and B2 spectra. Combinations with larger temperatures 
and lower transverse flow rapidities and vice versa lead to 
deuteron spectra with the same slope, but somewhat dif- 
ferent normalizations. The problem of different p and 
d slopes cannot be resolved in this way. Please note 
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that there is only a small temperature window for which 
the fugacity n/T needed to reproduce the normalization 
of the proton spectrum compares well with the value 
1.5 % [ijT % 2.0 from the experimental d/d ratio; due 
to the high nucleon mass, the nucleon and cluster yields 
are highly sensitive to the temperature. The fits tend to 
under- rather than overestimate the deuteron yields; this 
eliminates the simple suggestion that deuterons from the 
center of the fireball (which tend to have smaller Pj) are 
absorbed by rescattering, leaving only a thin surface shell 
contributing to the observed deuteron yield, because this 
would further aggravate the normalization problem. 

In a second step, we tested for fixed T = 140 Me V and 
/i/T = 1.75 different transverse flow proflles (2.7). For 
a = 0.5 and a — 2 good fits for the proton spectrum can 
be obtained with ry/ = 0.38 and r\f = 0.24, respectively. 
The effects of the modified flow proflle on the deuteron 
and _B2 spectra are small: lines no. 5 correspond to a = 

0. 5. which results in a weak, but insufficient rise of Bi 
with Mt; a — 1 (lines no. 7) made things worse. 

Finally, we fixed iJ./T=1.75 and tested transverse box 
profiles with different radii Pbox for the function H{R) at 
different lifetimes tq. As before, the remaining param- 
eters T and 77/ are adjusted to provide good fits to the 
proton spectrum. Lines no. 1 to no. 4 correspond to the 
following combinations of (t q, p box, T, rjf, fJ^/T), respec- 
tively (values of 77/ refer to (2.7) with fixed Ap — 7fm!): 

1. (7 fm/c, 10 fm, 142 MeV, 0.40, 1.75); 

2. (7 fm/c, 12 fm, 136 MeV, 0.34, 1.75); 

3. (8 fm/c, 12 fm, 134 MeV, 0.36, 1.75); and 

4. (9 fm/c, 14 fm, 126 MeV, 0.28, 1.75). 

All four data sets reproduce the slope of the preliminary 
deuteron and B2 data very nicely. They differ somewhat 
in the normalization of the deuteron spectrum. The pre- 
liminary data are located between curves 1 and 3, close 
to curve 2. The parameter choice no. 4 seems to be the 
most appropri ate w ith respect to the Gaussian parame- 

2(Ap) 



ter set of Sec. [IF, since pb 



14 fm for the 

box radius leads to the same transverse rms radius of the 
source (which is the observable determined from pion in- 
terferometry) as the Gaussian profile. However, it un- 
derestimates the yield of the deuterons and the value of 
B2, indicating too large an effective source volume for 
nucleons (see Sec. |IIE| and Eq. 3.2). Better results are 
obtained with somewhat smaller values of pbox and tq 
(curves 1-3 in Figs. ^,c). This is not ne cess arily in con- 
tradiction to the parameter set of Sec. 11 F, which was 
extracted from very central Pb-f-Pb collisions. For the 
less central NA44 collisions somewat smaller values of tq 
and/or Ap may indeed be expected. 

Given the preliminary nature of the data we did not 
make a big effort to optimize our fit to the data, e.g. 
make a combined fit to both proton and deuteron spec- 
tra. However, as minimum consistency requirements the 
lifetime tq and the flow velocity rjf must be large enough 
to allow for the transverse expansion from the cold Pb 
hard sphere radius (« 7fm) to the freeze-out radius pbox- 
This is an additional reason to reject parameter choice 



The outcome of this study is that Gaussian transverse 
density profiles cannot consistently reproduce both the 
proton and deuteron spectra, irrespective of the form of 
the transverse fiow profile. A transverse box profile for 
the density works quite well, both with linear (a = 1) 
and nonlinear (a = 0.5 and 2) transverse flow profiles. 
The source parameters needed to reproduce the deuteron 
yields and spectra are consistent with those extracted 
from HBT interferometry with pions. 

The physical interpretation of these results is that the 
data require more nucleons at higher transverse flow than 
the Gaussian density profile can provide even with steep 
velocity profiles. Other transverse density profiles like a 
Woods-Saxon or a doughnut profile may work similarly 
well or even better; on the other hand, the HBT ana- 
lysis of pion correlations has so far provided no hint for 
opaqueness of the source A good discussion of dif- 
ferent profiles, based on the analysis of E802 data from 
Si+Au collisions, can be found in |55|. 



C. Au+Au data from the AGS 

Although the size of a cold Au nucleus is about the 
same as that of a cold Pb nucleus, the lower beam energy 
at the AGS will req uire d ifferent fireball parameters than 
those given in Sec. H F , even if only collisions taken at 
the same 4% CTtot centrality are considered. A few simple 
checks with our theoretical results can nevertheless be 
performed. For E877, Johnson provides for the invariant 
coalescence factors in the lowest Pt bin the values = 
(1.5 ± 0.5) • 10-3 GeV2, Bs^ = (1.19 ± 0.29) • 10"*^ GeV^ 
and B3Hc = (1.25 ± 0.23) • 10"^^ GeV* Q. B^ is thus 
about a factor 2 larger than expected for central Pb-f Pb 
collisions at the SPS and seen by NA52. Using (6.2) and 
assuming that the correction factors (Ca) do not differ 
much between Pb-t-Pb and Au+Au collisions, we expect 



B2 (Au+Au) 
S2(Pb+Pb) 



Feff(Pb+Pb) 

Feff (Au+Au) 



P3 (Au+Au) 
S3(Pb+Pb) ■ 

(7.2) 



These rela tion s work out quite nicely when the Pb+Pb 
estimates (7.1) and the above experimental values for 
Au+Au are inserted. The E877 value of B2 is consistent 
with the preliminary results B2 = (1.8+1.0) • 10"^ GeV^ 
from E878 for central Au+Au collisions at y = 1.7 and 
B2 = (0.6 ± 0.4) • 10-3 GeV2 from E864 for 10% crtot 
Au+Pb coUisions at y = 1.9 p!]. 



VIII. SUMMARY 

We have presented an analytic treatment, supported 
by numerical checks, of the coalescence of two- and three- 
(anti)nucleon clusters in a relativistically expanding fire- 
ball. We showed that coalescence of nuclear clusters and 
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two-particle correlations between the final state momenta 
of pairs of identical particles can be characterized by the 
same set of effective source parameters, the "lengths of 
homogeneity". The only additional scale entering the 
coalescence probability is the intrinsic size of the clus- 
ter. We made the connection between coalescence and 
HBT interferometry explicit with the help of a simple, 
but phenomenologically successful parametrization of the 
expanding source which was already extensively studied 
in connection with two-pion HBT interferometry. For 
this model most calculations can be done in good ap- 
proximation analytically. 

Starting from the quantum mechanical definition of 
the deuteron momentum spectrum as a projection of 
the two-nucleon Wigner density of the fireball at nucleon 
freeze-out on the deuteron Wigner function, we used the 
formalism developed by Danielewicz and Schuck to 
bring the deuteron momentum spectrum into the form 
( |3.19 ) of a modified Cooper- Frye formula. Except for 
the "quantum mechanical correction factor" C{Rd,Pd) 
this formula looks exactly like a thermal spectrum of el- 
ementary particles with the mass of the deuteron, where 
the deuteron phase-space distribution is given by a prod- 
uct of thermal neutron and proton phase-space distribu- 
tions. The factor Cd is given by the integral (3.20) of 
the deuteron Wigner density over the homogeneity vol- 
ume of the nucleon source; for each deuteron momentum 
Pd it is peaked at the value Rd at which the flow ve- 
locity equals the deuteron velocity. We showed that the 
peak value for deuterons at rest in the fireball frame is a 
good approximation for (Cd) , the average of Cd over the 
freeze-out hypersurface, which is the relevant quantity 
for the shape of the deuteron momentum spectr um. ( Cd) 
can be expressed by the very simple formula ( [4.12| ) in 
terms of the ratios between the deuteron size d and the 
longitudinal and transverse lengths of homogeneity for 
nucleons at rest, 7?.||(m) and TZj^{m), respectively. In our 
expanding fireball model (C) turns out to be essentially 
momentum- independent for clusters of all sizes, and is of 
order (C) ~ 0.7 — 0.8 for 2- and 3-nucleon clusters formed 
in leOAGeV Pb-^Pb collisions. 

The shape of the deuteron moment um spectrum is thus 
given by the thermal model ansatz (2.3), with (C) lead- 
ing only to a modification of the normalization. For our 
Gaussian source parametrization the single particle spec- 
t rum c an be written, up to trivial factors, as a product 
( 2.14b| ) of an exponential i n M f with inverse slope ("ef- 
fective temparature" ) T* ( 2.16 ), a Gaussian in rapidity 
with width (Ar/)^/A, and an effective volume Vcs{A, Mt) 
whic h is again given by the homogeneity lengths , 
see ( ^.191 ). This allows to express the invariant coales- 
cence factor B2 fully in terms o f th e deuteron size and 
the lengths of homogeneity, see ( |6.3D . 

Unfortunately, the Gaussian model (Gaussian trans- 
verse density profile with linear transverse flow rapid- 
ity profile) leads to identical Aff-slopes for clusters of all 
sizes and therefore to Ba values which are essentially 
momentum- independent. This is contradicted by exper- 



iment. We found, in agreement with Polleri et al. [ p5[ , 
that the flatter deuteron than proton spectra require a 
density profile which gives for larger clusters more weight 
to regions of larger transverse flow. We found that a 
transverse box profile works very well and successfully re- 
produces the measured slope of preliminary data for 
the deuteron Mj-spectrum in Pb-|-Pb collisions. The con- 
sequence is a momentum-dependence of B2 which now 
rises as a function of Mt according to the simple gener- 
alized expression (6.4). 

It is interesting to observe that in this way the cluster 
spectra provide additional information about the source 
which cannot be extracted from HBT measurements. 
The latter constrain only the rms radii of the effective 
source, but not the shape of its spatial distribution. In 
the meantime, however, it has been found that a com- 
bined analysis of pion spectra and HBT correlations also 
prefers a transverse box profile over a Gaussian one be- 
cause the former gives a larger total pion yield than the 
latter, as required by the data [ [7l| . 
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TABLE I. The quantum mechanical correction factor C° 
for Hulthen and harmonic oscillator wave functions calculated 



with Eq. (3.21), for different fireball parameters at nucleon 
freeze-out (for details see text). 



TO [fm/c] 


9.0 


6.0 


T [MeV] 


168 


130 


100 


168 


130 


100 


Vf 


0.28 


0.35 


0.43 


0.28 


0.35 


0.43 


Hulthen 


0.86 


0.84 


0.80 


0.80 


0.78 


0.74 


harm. osc. 


0.84 


0.81 


0.76 


0.76 


0.72 


0.66 



19 



B,(d/p^) (GeV^) 

preliminary 

-2 



10 



-3 

10 



Pt=0, y=3,75 

p corrected for A decoy 



■k ■ 

J. : 

' our prediction 



B,(-d/p') (GeV^) ^ ^c-T (barn) 



-2 

10 



-3 

10 



Pt=0, y=3.75 

p corrected for A decoy 



, |.j 

|. f our prediction 



1 CTcuT (barn) 



FIG. 1. B2 for deuterons and antideuterons in 158AGeV Pb+Pb collisions as a function of the centrality, measured by 
NA52 Q. (The vertical scale has a systematic error of ~ 40%.) With al'^f fa 8.2 barn in Pb+Pb, the left-most data points 
correspond to a centrality of about 4%atot. The dashed lines show our prediction for B2 at 4%atot (see text). 
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FIG. 2. Fits to the proton mt-spectrum (a) and predicted deuteron Mt-spectra (b) and the ratio B2 = d/p^ (c) for different 
models for the transverse density profile. The analysis was performed with preliminary data (not shown) from the NA44 
Collaboration Q for 20%(Ttot Pb+Pb collisions at 158 AGeV in the rapidity range 1.9 < y < 2.3. Lines 1-4: box profile for the 
transverse density with linear transverse flow proflle; lines 5-7: Gaussian transverse density profile with power law transverse 
flow profiles, with exponents a = 0.5, a = 1.0 and a = 2.0. The preliminary data lie close to line 2. For details see text. 
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